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ABSTRACT 

Context. The nuclear star cluster of the Galaxy is an important template for understanding its extragalactic counterparts, which can 
currently not be resolved into individual stars. Important drawbacks of observations of the Galactic center are, however, the presence of 
strong and spatially highly variable interstellar extinction and extreme crowding of the sources, which makes the use of adaptive optics 
techniques necessary. Both points pose serious obstacles to precise photometry that is needed for analyzing the stellar population. 
Aims. The aims of this work are to provide accurate photometry in multiple near-infrared broadband filters, to determine the power- 
law index of the extinction-law toward the central parsec of the Galaxy, to provide measurements of the absolute extinction toward 
the Galactic center, and finally to measure the spatial variability of extinction on arcsecond scales. 

Methods. We use observations of the central parsec of the Milky Way that were obtained with the near-infrared camera and adaptive 
optics system NAOS/CONICA at the ESO VLT unit telescope 4. The photometric method takes into account anisoplanatic effects and 
limits the corresponding systematic uncertainties to < 2%. Absolute values for the extinction in the H, Ks, and L' -bands as well as 
of the power-law indices of the H to Ks and Ks to L' extinction-laws are measured based on the well-known properties of red clump 
stars. Extinction maps are derived based on H - Ks and Ks - L' colors. 

Results. We present Ks-band photometry for ~ 7700 stars, and additionally photometry for stars detected in the H and/or L' -bands. 
From a number of recently published values we compute a mean distance of the Galactic center of R = 8.03 ± 0.15 kpc, which has 
an uncertainty of just 2%. Based on this Rq and on the RC method, we derive absolute mean extinction values toward the central 
parsec of the Galaxy of A H = 4.48 ± 0.13mag, A Ks = 2.54 ± 0.12mag, and A L i = 1.27 ± 0.18mag. We estimate values of the 
power-law indices of the extinction-law of a H - Ks = 2.21 ± 0.24 and a Ks -u = 1.34 ± 0.29. A A"i-band extinction map for the Galactic 
center is computed based on this extinction law and on stellar H - Ks colors. Both its statistical and systematic uncertainties are 
estimated to be < 10%. Extinction in this map derived from stellar color excesses is found to vary on arcsecond scales, with a 
mean value of A Ks = 2.1 A ± 0.30 mag. Mean extinction values in a circular region with 0.5" radius centered on Sagittarius A* are 
Ans g rA, = 4.35 ± 0.12, A Ks j g ,-A, = 2.46 ± 0.03, and A UtSgrA , = 1.23 ± 0.08. 



1. Introduction 

Nuclear star clusters (NSCs) are located at the photometric and 
dynamical centers of almost all spiral galaxies. Their detailed 
study has only become possible via HST observations in the 
1990s, but they are of great interest because they are the densest 
and most massive clusters in their host galaxies. Frequently they 
contain supermassive black holes at their centers. They show 
signs of recurring star formation, and the discovery of rotation of 
the NSCs of NGC4244 (|Sefh et al.|2008|> and of the Milky Way 
([Trippe et aTJ|2"0"08j |Schodel et al'||2009|l indicates that NSCs 



may at least partially form by accretion of material from their 
host galaxy. A concise review of the properties of NSCs and hy- 
potheses on their formation is given in Boker (20081). 

Studying the NSC at the Galactic center (GC) is of great im- 
portance for the understanding of NSCs in general because in 



the next decades it will remain impossible to resolve any but 
the brightest individual stars in extragalactic NSCs, even with 
50m-class telescopes. A further reason for the importance of the 
GC NSC is that the existence of the supermassive black hole, 
Sagittarius A*, at the GC is well established and the mass of this 



black hole has been determined with high precision (e.g., Ghez 
|et al.||2008[ |Gillessen et al.||2009| >. Strong interest in the stellar 
population at the GC was raised by spectroscopic observations, 
which have revealed the existence of numerous massive, young 
stars within 0.5 pc of Sagittarius A*. Those young stars are wit- 
nesses of a starburst event that occurred about 4 million years 
ago. The young stars are arranged in a non-isotropic pattern and 
at least 50% of them are located in a rotating disk (see |Bartko| 
|etal.|2009t|Luet al.|2009] >. 
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Unfortunately, there are two important drawbacks in obser- 
vational studies of the Milky Way's NSC: extreme crowding and 
high extinction. While the great progress made in adaptive optics 
techniques at 8-10m-class telescopes in the past years has been 
highly successful in reducing the problem of crowding, it nev- 
ertheless remains a serious issue. Even with instruments such as 
NACO at the ESO VLT, strong crowding limits the point-source 
sensitivity within ~ 0.5 pc of Sagittarius A* so significantly that 
only stars more massive than ~ 2 M can be detected in this re- 
gion. This means that we can still detect only of the order 10% of 
the stars in the NSC ( |Schodel et al.|2007a| l. Additionally, it is no- 
toriously difficult to extract accurate photometric measurements 
from AO observations of crowded fields, particularly when the 
field-of-view exceeds the size of the isoplanatic angle. 

As concerns extinction, we can consider ourselves lucky that 
we can observe the GC at all. Within 2-3 pc of Sagittarius A* 
there exist plenty of clouds that are almost or completely opaque 
even in the near-infrared (for example, see Fig. 1 in Schodel 
|et al. |2007a[ ). Average extinction toward the central parsec in 
the K-band reaches Ak ~ 3. As concerns extinction at vi- 

9, 
16 



sual wavelengths, |Rieke & Lebofsky 



concerns 
1985) found A V IA K 



but more recent studies indicate higher values o f A V /A K 
(|Nishiyama et al.| |2008 ), or even A V /A K ~ 29 ( |Gosling eFaT 
2009 ). Therefore, the stellar population at the GC becomes ob- 



servable only at near-infrared wavelengths. Only the brightest 
stars can be seen in the /-band. Reasonable S/N on the stars 
is only reached in the H, K, and L bands. This restriction to 
a tight wavelength window, where most stars show only small 
color indices, combined with the fact that the extinction toward 
the GC is not only high, but also variable on small spatial scales 
of a few arcseconds ( |Scoville et al.|[2003| |Gosling et al.| [2006 ; 
Schodel et al. |2007a I, makes photometric studies of the stellar 



population at the GC extremely difficult 

Although adaptive optics assisted integral field spectroscopy 
has resulted in the classification of hundreds of stars within 
~ 0.5 pc of Sagittarius A* (e.g., IManess et al.||2007[ IPaumard 



|et al.||2006t |Do et a l. 2009a), this technique is extremely time 
consuming in both observations and data analysis. In order to 
efficiently probe large fields in the Galactic center for the pres- 
ence of young, massive stars, it is therefore desirable to improve 
the photometric observations in a way that enables us to at least 
crudely distinguish whether a star in the GC is of early or of late 
spectral type. How this can be done in spite of the difficulties de- 



scribed above has been demonstrated recently by Buchholz et al. 



(2009). They used intermediate band filters across the K-band, 
calibrated photometry locally with the help of red clump stars, 
and used the CO bandhead absorption feature of late-type stars 
as a distinguishing criterion. 

Accurate knowledge of extinction has impact in many dif- 
ferent aspects of GC research: the stellar surface number den- 
sity is not only a function of distance from Sagittarius A*, but is 
also dependent on local extinction ( Sch odel et al.|2007a) >; extinc- 
tion is an important ingredient for modeling stellar atmospheres 
(e.g., Najarro et al.|1 997; Martins et al. 2007); extinction in the 



infrared is also important for estimating the intrinsic brightness 
of the near-infrared emission from the accretion flow or outflow 
related to Sagittarius A* (e.g., |Genzel et al. 2003; Eckart et al. 
|2006|[Do et al.|2009b[ ). The mentioned examples show that one 
must know not only the absolute value of extinction, but also 
how it changes as a function of wavelength and how it varies on 
spatial scales as small as one to a few arcseconds. 

Until just about a decade ago, the extinction law from near- 
to mid-infrared wavelengths (~ 1 — 8 /mi) was thought to be 
described well by a "universal" power-law, Aj oc A~ a , with a 



power-law index of a x 1.75 (e.g., |Draine||1989] ). First doubts 
were shed on this model by spectroscopic observations of the 
GC with ISO SWS that showed clear deviations from this model 
at wavelengths > 3 jum ( |Lutz et ai][l9 96 ). Subsequent studies 
have confirmed that the near- to mid-infrared extinction-law is 
more complex than a simple power-law. While it seems well 
established by now that a simple power-law is an adequate de- 
scription of the wavelength dependence of extinction in the near- 
infrared, albeit with a stronger wavelength dependence as previ- 
ously thought, the situation is less clear longward of ~ 3 /mi. 
Recent investigation shows that the extinction law decreases 
much less rapidly at A > 3 /mi than suggestd by earlier estimates 
and may actually be almost constant at 5 - 10 /mi (see Nishiyama 
et al.|200"9] for an extensive discussion of these points). Taking 
the problem even a step further, in a recent study of extinction 
toward the nuclear bulge Gosling et al. (2009 ) find that the index 
of the power-law extinction-law can change on spatial scales as 
small as 5 arcseconds. This signifies that we may be forced to 
abandon the concept of a universal extinction-law. 

It is evident that the topic of extinction toward the central 
parsec of the Galaxy needs to be re-visited. Extinction maps for 
the central 0.5-1 pc were recently presented by Sch odel et al.| 
(2007a) and Buchholz et al. (2009 ). However, these papers were 
not mainly focused on extinction and the extinction maps were 
rather side-products. Both papers used the old Aj oc A~ lJ5 
extinction-law ( |Draine| [l989). Neither the extinction law nor 
absolute extinction were examined in depth. Also, the abso- 
lute photometric calibration in our previous works was of lim- 
ited accuracy because the PSF for the stars was either truncated 
( Schodel et al.|20O7a| l or because complex local calibration pro- 
cedures had to be adopted ( jBuchholz et al.|2009) >. 

The aim of this paper is to carry us a step forward in obser- 
vational studies of the GC NSC by addressing two of the above 
described problems: accurate photometry with AO over a large 
FOV and accurate measurements of extinction. Schoedel (A&A, 
2009, accepted for publication) has recently presented a method 
for the photometric analysis of AO data with spatially variable 
PSF, but a very limited number of suitable PSF reference stars. 
Based on this method, in this paper we analyze imaging obser- 
vations of the GC in the H, Ks, and L' bands and provide ac- 
curate photometry (with both statistical and systematic uncer- 
tainties of the order ~ 5% at Ks) for several thousand sources 
observed in these bands within ~1 pc of Sagittarius A*. We de- 
termine both the absolute extinction and the extinction law at the 
various bands and provide a detailed extinction map (resolution 
of 1" - 2") for the central parsec of the GC. 



2. Observations 

The observations used in this work were obtained with 
the near-infrared camera and adaptive optics (AO) system 
NAOS/CONICA (short NACO) at the ESO VLT unit telescope 4 
Q The Ks ss 6.5 - 7.0 supergiant IRS 7 was used to close the 
loop of the AO, using the NIR wavefront sensor. The sky back- 
ground was measured on a largely empty patch of sky, a dark 
cloud about 400" north and 713" east of the target. Sky subtrac- 
tion, bad pixel correction, and flat fielding were applied to the 
individual exposures. The NACO S27 camera, with a pixel scale 
of 0.02705" /pix, was used for both H and Ks band observations. 
The field-of-view (FOV) of a single exposure is thus 28" x 28". 



1 Based on observations made with ESO Telescopes at the La Silla 
or Paranal Observatories under programmes ID 077.B-0014 and 081. B- 
0648 
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Fig. 1. Mosaic image of the Ks-band observations. North is 
up and east is to the left. Offsets are given in arcseconds from 
SgrA*. The guide star, IRS 7, is the brightest star within the 
rectangle, about 5.5" north of SgrA*. The black rectangle in- 
dicates the area of overlap between the 4 dither positions. Blue 
crosses mark early-type stars identified by their H - Ks colors 
with extinction/reddening correction applied (see section|5] and 
Fig.©. 
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Fig. 2. Mosaic image of the L'-band observations. North is up 
and east is to the left. Offsets are given in arcseconds from 
SgrA*. 



The observations were dithered by applying a rectangular dither 
pattern with the center of the dithered exposures positioned 
approximately at (8.0", -2.6"), (-6.1", -2.7"), (-6.1", 11.2"), 



and (8.1", 11.3") east and north of SgrA*. In the text we refer 
to these four offsets as dither positions 1, 2, 3, and 4. The com- 
bined FOV of the observations is about 40" x 40" and is offset 
to the north with respect to Sgr A* because the guide star IRS 7 
is located about 5.6" north of SgrA*. 

Seeing in the optical ranged between 0.6" and 1.0" for 
both H and Ks observations. According to the ESO Ambient 
Conditions Database^jsky transparency variations during the ob- 
servations were below 1%. The achieved Strehl ratio ranged be- 
tween ~ 20% (~ 15%) near the guide star and ~ 8% (~ 8%) 
at 25" distance from the guide star in the Ks-{H-)ban&. The 
Strehl ratio was estimated using the Strehl algorithm of the ESO 
eclipse software package (Devilla rd|1997 1 on PSFs extracted at 
various positions in the image. From the multiple measurements 
we estimate the 1 <x uncertainty of the measured Strehl ratio to 
~ 3%. TableQ] summarizes the observations. The detector inte- 
gration time (DIT) was set to 2.0 s in order to avoid saturation of 
the brightest stars. The central pixels of a small number of the 
brightest stars were saturated in Ks. Their cores were repaired 
using nearby unsaturated stars and the PSF extraction routine of 
StarFinder. After 28 DITs, the instrument averaged the data to 
a single exposure (NDIT= 28). In this way, 8 individual expo- 
sures were obtained per dither position. The exposures of each 
respective dither position were aligned (to compensate for small 
residual shifts) with the jitter algorithm of the ESO eclipse soft- 
ware package (Devillard 1997). We show the combined FOV 
of the Ks-band observations in Fig.[T] Note that the photometry 
and astrometry in this work was done on the combined images 
of each dither position and not on the combined mosaic of all 
images (as shown in Fig.[T]i in order to have a constant signal-to- 
noise ratio over the entire images. The ~ 13.5" x 13.5" overlap 
area between the four dither positions is indicated by the rectan- 
gle in Fig.[T] 

Seeing in the optical ranged between 1.0" and 2.0" for the 
U observations. The NACO L27 camera, with a pixel scale of 
0. 02719" /pix, was used for the L' observations. According to 
the ESO Ambient Conditions Database sky transparency varia- 
tions during the observations were below 1%. The observations 
were random jittered, with a jitter box width of 20". Offsets by 
60" in random directions were alternated with observations of 
the target. The frames taken on the offsets frames, where the 
stellar density is lower, were used to construct sky frames. The 
7 offset frames nearest in time to a given exposure on the tar- 
get were median combined and subtracted as sky measurement. 
Flat fielding and bad pixel correction were applied. Finally, all 
exposures were combined with the jitter algorithm of the ESO 
eclipse software package. The resulting mosaic image is shown 
in Fig.|2] 



3. Photometry 

Due to the extreme source density at the GC aperture photome- 
try cannot be applied accurately to GC imaging data. Therefore 
PSF fitting techniques must be used. In AO imaging anisopla- 
natic effects can lead to a systematic bias in the photometry of 
point sources when PSF fitting is used for estimating the flux of 
the stars. This bias depends on various factors such as distance 
from the guide star, observing wavelength, and atmospheric con- 
ditions at the time of the observations. In case of the //-band 
data used in this work the systematic error of the photometry 
can reach 0.25 mag across the FOV, if a single, non-variable 



http://archive.eso.org/asm/ambient-server?site=paranal 
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Table 1. Details of the imaging observations used in this work. 



Date 


■^central 


AA 


N" 


NDIT" 


DIT C 






[jum] 






[s] 


29 April 2006 


1.66 


0.33 


32 


28 


2 


29 April 2006 


2.18 


0.35 


32 


28 


2 


26 May 2008 


3.80 


0.62 


202 


150 


0.2 



a Number of (dithered) exposures 

b Number of integrations that were averaged on-line by the read-out 
electronics 

c Detector integration time. The total integration time of each obser- 
vation amounts to NxNDITxDIT. 



PSF is used (see Schoedel, A&A, 2009, accepted for publica- 
tion). For example, anisoplanatic effects make stars at distances 
larger than the isoplanatic angle appear elongated, with their 
major axes pointing toward the guide star. This can be nicely 
seen in the image shown in Fig. [2] (note: usually anisoplanatic 
effects are stronger at shorter wavelengths, but in the observa- 
tions used for this paper, it happens that those effects are stronger 
in the L'-band than in the Ks-band image). Various methods 
for dealing with anisoplanatic effects have been suggested. A 
straightforward and effective method to deal with anisoplanatic 
effects is presented in a companion paper (Schoedel, A&A, 
2009, accepted for publication). It not only provides accurate 
point source photometry but also reliable estimation of the dif- 
fuse background emission. In brief, the method consists of two 
steps: (a) Wiener deconvolution, using a point spread function 
(PSF) extracted from the guide star (or bright stars near the 
guide star), in this case IRS 7; (b) photometry and astrometry 
on the deconvolved image via PSF fitting with the StarFinder 
( |Diolaiti et al.|2 000) program package. The image is partitioned 
into overlapping sub-frames that are smaller than the isoplanatic 
angle. Those sub-frames are analyzed independently. As shown 
in Schoedel (A&A, 2009, accepted for publication) this method 
minimizes the systematic uncertainty of the point source pho- 
tometry due to anisoplanatic effects to < 0.02 mag across the 
28" x 28" FOV of the images at each dither position. 

After deconvolution with PSFs extracted from IRS 7, the H- 
and Ks-band images from the four dither positions were parti- 
tioned into 13x13 overlapping sub-frames of 256 x 256 pixel 2 
(~ 6.9" x 6.9"). The shifts between sub-frames were thus just 
64 pixels in x- and/or y-direction to assure large overlap between 
the sub-frames. We used the following StarFinder parameters: 
minjoorr — 0.9, thresh — [5., 5.], and backJbox — 30. In or- 
der to avoid spurious detections, point source detection was also 
run on non-deconvolved images and only sources that were dis- 
covered in both the deconvolved and non-deconvolved images 
were included in the final list of detected stars. The images for 
the diffuse background and residuals were obtained by recom- 
bining the corresponding sub-frames. The method provides two 
kinds of uncertainties for the photometric and astrometric mea- 
surements: formal uncertainties, computed by the StarFinder al- 
gorithm from Gaussian and photon noise, and PSF uncertain- 
ties, resulting from incomplete knowledge of the local PSF. The 
latter is calculated by comparing multiple measurements of the 
stars because they are present in several sub-frames. For stars 
without multiple measurements we adopted a PSF uncertainty 
of 0.02 mag for the //-band and 0.015 mag for the Ks-band 
(see Fig. [3]). The results of multiple measurements were aver- 
aged. Since deconvolution changes the noise properties of the 
image, the formal uncertainties computed by StarFinder are too 
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Fig. 3. Statistical photometric uncertainty vs. magnitude for the 
H- (top), Ks- (middle), and L '- (bottom) bands. 



small (for details see Schoedel,A&A, 2009, accepted for publi- 
cation). Monte Carlo simulations were used to determine a cor- 
rection factor with which we scaled the formal uncertainties. The 
final uncertainties were calculated by quadratically combining 
the re-scaled formal fit uncertainties with the PSF uncertainties. 
The data in the overlap region between the four dither positions 
were averaged. Stars present in the overlap region between the 
dithered fields in the H and Ks observations were used to check 
the photometric accuracy. Photometry on the stars in the over- 
lapping fields shows that the uncertainties estimated by the ap- 
plied method are correct (see Schoedel, A&A, 2009, accepted 
for publication). 
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Fig. 4. Color-magnitude diagrams, with Ks vs. H - Ks shown 
in the top panel and L' vs. Ks - L' in the middle one. Bottom: 
Completeness (due to crowding) corrected L', Ks, and H lumi- 
nosity functions. The dashed green lines are fits with an expo- 
nential plus a Gaussian. 



Photometry on the L' image was performed in a very similar 
way, except that there was only one combined frame to be ex- 
amined. The PSF for initial linear deconvolution was extracted 
from IRS 7. The size of the sub-frames for local PSF fitting was 



300 x 300 pixel 2 , with a shift of 50 pixels between the frames. 
Since the L' data are very rich in structure as concerns the diffuse 
emission, we set the backJbox to a smaller value (20 pixels). 

Zero points for the NACO instrument for various combina- 
tions of cameras, filters, and dichroics (that split the light be- 
tween wavefront sensor and camera) are determined routinely 
within the ESO instrument calibration plan. Zero points for all 
filters used in this work and for the corresponding setup (cam- 
era S27, dichroic N20C80) were determined via observation of 
a standard star during the same night as the observations: ZPh = 
23.64 + 0.05, ZP Ks = 22.85 + 0.05, and ZP L , = 22.38 + 0.15. 
Additional uncertainty of the zero point is caused by the uncer- 
tainty of the PSF that is used for deconvolution. From experi- 
ments (choosing different PSF sizes) we estimate this additional 
uncertainty to < 0.03 mag. We combine this additional uncer- 
tainty quadratically with the previous uncertainties and obtain 
final zero point uncertainties of AZP# = 0.06, AZPk s = 0.06, 
and hZPu = 0.15. The annual mean and standard deviation of 
the zero points corresponding to the year and the setups used for 
the observations in this work are ZPn,mean,im> — 23.61 ±0.10, 
ZP Ks sneanam = 22.82 + 0.07, and ZP z ,- jmean(20 08 = 22.25 ± 0.17. 
This agrees well with the zero points and corresponding 1 cr used 
for the photometry in this work. 

Plots of photometric uncertainty vs. magnitude for all point 
sources discovered in the individual bands are shown in Fig. [3] 
We combined the photometric results to a final list, requiring that 
each star must be detected in at least two different bands. The 
final lists of detected point sources was aligned with the coordi- 
nate system of the /^-observations. The full list of all stars and 



their fluxes in the different bands is given in Table A. 2 available 
online. 

Because of the non-standard photometric methodology ap- 
plied here (deconvolution and extraction of local PSFs) it was 
not straightforward to determine incompleteness due to crowd- 
ing via the standard method of inserting and recovering artif- 
ical stars. The main difficulty here is the change of the PSF 
across the FOV. We therefore estimated completeness directly 
from the data using the methodolog y suggested by |Eisenhauer| 
et al. ( 1998 1 and further described in Harayama et al.| (2008 ). In 
brief, it consists of deriving a critical distance for each possible 
magnitude difference between two given stars. A star closer to a 
brighter one than the critical distance will not be detected. Once 
the critical distance is known for all relevant magnitude differ- 
ences, completeness maps can be created. Table|2]gives the com- 
pleteness of our data for different magnitudes and bands. Some 



of the assumptions of the methodology suggested by |Eisenhauer 
et al. ( | 1998 1 may be violated in our case, particularly becase 



the stellar density in the GC FOV changes strongly because of 
effects of extinction and cluster structure. For this reason we es- 
timated the critical distances in various ways, e.g. by analysing 
different parts of the FOV. In Table|2]we therefore give two val- 
ues for estimated completeness, a high and a low one, for each 
magnitude and band. Fortunately it turns out (see following sec- 
tion) that completeness correction does not have any significant 
influence on our results. 



4. Constraining extinction with the red clump 

Color-magnitude diagrams of our data are shown in Fig.|4 The 
densely populated area around Ks as 15.6 and L' « 14.2 is 
due to red clump (RC) stars, which are the equivalent to hori- 
zontal branch stars in metal-rich populations. The luminosities 
and colors of RC stars are narrowly distributed and are only 
weakly dependent on age and metallicity, particularly for clus- 
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Table 2. Completeness levels in percent for given bands and 
magnitudes. 



mag 


H 


Ks 


L' 


12 


100-100" 


98-100" 


93-96" 


13 


100-100 


96-98 


84-91 


14 


99-99 


91-97 


64-77 


15 


97-99 


83-92 


54-69 


16 


93-97 


59-82 


42-58 


17 


83-91 


47-70 


n.a. 


18 


66-82 


28-50 


n.a. 


19 


50-70 


n.a. 


n.a. 



A high and a low value of estimated completeness are given for 
each combination of band and magnitude. 



ters older than -1.6 Gyr (e.g.,|Castellani et al.|1992||Alves|2000 
Grocholsk i & Sarajedini 2002 1. 

The completeness corrected luminosity functions (LFs) in 
the H, Ks, and L' bands are shown in the bottom panel of Fig. [4] 
The RC feature is clearly visible as a bump in all LFs. We fitted 
all LFs with the combination of exponentials plus Gaussians. 
The locations of the peaks of the RC as derived from these 
fits are H RC = 17.60 ± 0.04(17.58 + 0.016), Ks RC = 15.59 ± 
0.04(15.57 ± 0.009), and L' RC = 14.25 ± 0.03 (14.22 ± 0.008). 
Values in brackets are the results obtained if no crowding com- 
pleteness correction is applied. The uncertainties of the fits are 
smaller than the uncertainty of the absolute photometric calibra- 
tion (0.06 mag for H and Ks, and 0.15 mag for L'), In the follow- 
ing analysis we will use the numbers based on the completeness 
corrected LFs. We note however, that this crowding complete- 
ness corection has a negligible effect, of the order 10% of the 
uncertainty, on the final results. We attribute this to two reasons: 
(a) the RC is located in a magnitude range where completeness is 
fairly high; (b) the dominating uncertainty is the one of the abso- 
lute photometric calibration. That incompleteness due to crowd- 
ing is not a serious issue can also be seen when comparing the 
KsLF from this paper with the one shown in Fig. 10 of Schodel 
et al.|p007a[ ). A different, slightly deeper data set was used for 
the latter work, but we can see that the RC bump in the KsLF is 
negligibly affected by incompleteness. Note, however, that abso- 
lute calibration was not a goal of Schodel et al. ( 2007a i, there- 
fore, there RC peak is shifted by ~ 0.3 mag to brighter magni- 
tudes. Finally, we may ask whether the strong extinction toward 
the GC may affect the RC peak in the LFs of the different fil- 
ters and thus bias the peak toward brighter values. RC stars are, 
however, well within the sensitivity of the observing setup used, 
at all bands. For example, an isolated Ks = 20 source would 
be detected with NACO and the settings used for this work at a 
level of ~ 10 o~. The only incompleteness that affects us here is 
the one due to crowding, which has fortunately a negligble effect 
on our results. 

Since RC stars can be found in large numbers toward the 
Galactic bulge and center, they are highly useful for studying 
Galactic extinction and structure (e.g., Wozniak & Stanek 1996; 
Udalski|2003||Sumi|2004]|Nishiyama et al.|2005||2006a|b| l. The 



absolute magnitude of RC stars has been estimated to M K = 
-1.54 ± 0.04 (in the 2MASS s ystem) by |Groenewegeri| p508] > 



and to M K = -1.61 ± 0.03 by |Alves| pOOO) . Both values are 



based on the Hipparcos catalogue. Here we adopt the more 
recent value for M K - The difference between the K and Ks- 



magnitude is » 0.01 mag (see discussion in Nishiyama et al. 
2006b|. [Nishiyama et al.| ( |2006b| estimated the distance to the 




2 2.2 
mag H -mag Ks 

Fig. 5. Peak K s-magnitude vs. H - Ks color of the red clump. 



GC based on photometry of RC stars, calculating the distance 
modulus as 



(m - M) = Ks RC ; mtI . - M Ks + AM K , 



(1) 



where ^^Rc,intr. is the intrinsic, reddening free magnitude of RC 
stars, Mks is their absolute /Ts-magnitude, and AMk the popula- 
tion correction for M Ks . We adopt the same values as Nishiyama 
[eTaTI ( |2006b"l ) for AM K = -0.07 and use M Ks = -1.54 (see 
above). By assuming a GC distance, we can use the above equa- 
tion to estimate the reddening free magnitude of RC stars. By 
comparing it subsequently with the magnitude of the RC clump 
peak as inferred from the luminosity function, we can then ob- 
tain an estimate of the absolute extinction in Ks toward the GC. 

Numerous efforts to determine the distance of the GC have 
been made recently (the cited uncertainties are the quadrati- 
cally combined statistical and systematic uncertainties given in 
the respective works): 7.52 ± 0.36 kpc (bulge red clump stars, 
2006b|l, 8.33 + 0.3 5 kpc (stellar orbit around 

0.6 kpc (stellar or- 
0.45 kpc 



Nishiyama et al 



Sagittarius A*, |Gillessen et al.||2009), 8.0 : 
bit around Sagittarius A*, |Ghez et aLl[2008)l, 7.94 



(Population II Cepheids and RR Lyrae stars Groenewegen et al. 
2008), 8.07 + 0.35 kpc (proper motion and line-of-sight velocity 
dispersion of late-type stars in central parsec, Trippe et al. 2008), 
8.4 ± 0.6 kpc (trigonometric parallaxes of masers measured with 
VLBI,|Reid et al.|2009a| , 8.24 + 0.43 kpc (Mira stars, |Matsunaga| 
|et al.|2009| >, and 7.9 + 0.8 kpc (trigonometric parallax of Sgr B2 
measured with VLBI, |Reid et aT1|2009b| . These distance esti- 
mates are based on different data sets and different methods and 
can probably be considered largely independent measurements. 
It appears therefore justified to combine these measurements to 
a weighted mean of Rq = 8.03 ± 0.15 kpc. This corresponds to a 
distance modulus of 14.52 + 0.04. This simple exercise demon- 
strates that - under the assumption that these measurements are 
indeed statistically independent - Rq is already known with the 
surprisingly small uncertainty of just 2%. Using this value we 
obtain a reddening free magnitude of the RC bump at the GC of 
^•SRQintr. = 13.05 + 0.09. All uncertainties have been quadrat- 
ically combined (A(m - M)o = 0.04 mag, AM Ks = 0.04 mag, 
AAM K = -0.07 mag). 

The fit of the RC bump in the KsLF with a Gaussian (plus 
an exponential to fit the underlying luminosity function) results 
in Ks R c = 15.59 + 0.07. Here, the uncertainties from the fit of 
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Fig. 6. Top left: Map of A Ks for the GC derived from H - Ks colors of all stars detected in these bands, after exclusion of fore- and 
background stars. Top right: Map of the resolution of the extinction, i.e. the radius in arcseconds at each position within which the 
20 stars for the given measurement were detected. Bottom left: Map of Ak s for the GC derived from Ks - L' colors. Bottom right: 
Map of the statistical uncertainty of A Ks , corresponding to the extinction map shown top left. 



the RC peak and from the absolute photometric calibration have 
been combined quadratically. From the difference with the theo- 
retical reddening free magnitude of the RC bump, this results in 
an average extinction of Ak s = 2.54 ± 0.12 toward the GC. 

Assuming H - Ks = 0.07 ± 0.03 and Ks - U = 0.07 ± 
0.03 for the RC stars, we obtain in the same way estimates for 
A H = 4.48 ±0.13 and Au = 1.27 ± 0.18. Note that these values 
depend only on the well-known properties of RC stars and on 
the GC distance. They are independent of any assumptions on 
the extinction-law. 

Assuming the validity of a power-law for the extinction law 
between H and Ks and between Ks and L', respectively (see, 
e.g.,|Nishiyam a et al.|2009| l, we can derive from the above mea- 
surements the respective power-law indices Qh-ks = 2.21 ±0.24 



and oks-u = 1-34 ± 0.29. Here the uncertainties include the 
uncertainties of the effective wavelengths at the correspond- 
ing bands. The effective wavelengths adopted in this work are 
^ ff , H = 1.677 ± 0.018 //m, i effjKs = 2.168 ± 0.012/an, and 
'leff.L' = 3.636 ± 0.01 2 fim. Details on how these effective wave- 
lengths have been calculated are given in appendix[A| 

As can be seen in the color-magnitude diagram shown in 
Fig.|4](top left) there is a clear dependency of the Ks peak mag- 
nitude of the RC feature on H - Ks. Because of the narrow in- 
trinsic distribution of colors and magnitudes of the RC stars we 
expect this trend to be almost exclusively due to differential ex- 
tinction. In order to quantify this relation, we extracted the sLF 
in narrow ranges of H - Ks. The peaks of the RC features were 
then fitted with a Gaussian. This method is very similar to the 
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one applied by Nishiyama et al. ( 2006a ), who examined the loca- 
tion of the RC feature in different fields toward the GC in order to 
derive the reddening law. Figure[5] shows the thus measured RC 
Ks peak magnitude vs. H-Ks along with a linear fit. Assuming 
A A oc A~ a "- Ks , the slope of the line gives au-Ks = 2.54 + 0.15. 
This value agrees within the uncertainties with the one derived 
from the RC bump in the H and Ks luminosity functions (see 
above) and thus serves as a cross-check. However, we consider 
the value au-Ks = 2.21 ± 0.24 derived from the RC as more 
reliable (see discussion section). 

The Ks vs. H-Ks color magnitude diagram shows a decreas- 
ing completeness with increasing reddening. This is reflected in 
Fig. [4] by the increasing faint magnitude limit of the cloud of 
dots toward redder colors. We have seen above that incomplete- 
ness does not significantly affect the measured peak of the RC 
in the KsLF. We have not used any completeness correction for 
the fit of the RC peak in the color magnitude diagram, but tested 
its reliability by omitting the two "reddest", i.e. highest H - Ks, 
data points. The fit then results in cvh-ks = 2.71 +0.18. If we ad- 
ditionally omit the bluest data point from the fit, which may be 
an outlier, we obtain uh-Ks = 2.45 ± 0.14. We conclude that in- 
completeness due to reddening does not have any significant in- 
fluence on the value of (Xh-ks derived from the color-magnitude 
diagram in Fig. |5] 

As can be seen in the upper right panel of Fig. [4] there is 
no clear trend of the RC visible in the plot of L' vs. Ks - U . 
This is mainly because the number of stars measured in L' is 
smaller (smaller FOV) and because the range spanned by the 
Ks—L' colors is limited due to the weak extinction in U . We can 
therefore not apply this method to determine a^s-L' (see below 
for a different method, however). 



5. An extinction map for the GC 

We produced an extinction map for the FOV of the H and Ks ob- 
servations of the GC by using the median of the H-Ks colors of 
the 20 nearest stars at each position. For this purpose we used all 
stars from our list, but excluded extremely blue (H - Ks < 1 .8) 
or red (H - Ks > 2.8) ones as fore- or background stars. The in- 
trinsic colors of almost all stellar types observable at the GC are 
similar at these wavelengths. Nevertheless, we tried to improve 
the estimates by applying a 0th order correction for different in- 
trinsic stellar colors. Stars fainter than Ks = 14.5 were assigned 
(H-Ks)q = 0.07 and stars brighter than this value were assigned 
(H - Ks)o = 0.2. Known early- type stars (taken from Buchholz 
et al. 2009[ ) were assigned (H - Ks)q = -0.03. In order to assign 
these intrinsic colors, we used the GC distance, approximate ex- 



be crowded. For example, the majority of stars within 0.5'' 
Sgr A* are of early-type (e.g. Eisenhauer et al.|2 005 ) 



of 



tinction, and the magnitude of the stars to roughly guess their 



type (see, e.g., Table 3 in Buchholz et al. 


2009) and assigned 


colors from tables 7.6 and 7.7 in|Cox||2000 





In order to check the importance of distinguishing different 
intrinsic stellar colors we also created a map of H - Ks assum- 
ing just one uniform intrinsic color, H — Ks — 0.07 for all stars. 
The resulting map looks almost identical to the one including the 
more refined corrections (the mean difference between the two 
maps is A(H - Ks) - 0.02 ± 0.02). We conclude that it would 
be, in principle, acceptable to assume just one common intrin- 
sic color for our entire sample of stars because only ~ 16% of 
our sample of stars are brighter than Ks = 14.5. The majority 
are, in fact, red clump stars. Although this shows that distin- 
guishing different types of stars in the 0th order color correction 
is almost negligible, we nevertheless include this correction in 
order to avoid bias in regions where certain types of stars may 



The map of color excesses, E(H - Ks), was converted to 
an extinction map by assuming a power-law extinction-law with 
an index an-Ks = 2.21 ± 0.24 for the effective wavelengths of 
A H = 1.677 ±0.018//m and A Ks = 2.168 ±0.012//m. The uncer- 
tainty of these three parameters combined with the uncertainty 
of the photometric zero points and of the intrinsic stellar colors 
results in a combined systematic uncertainty of ~ 9% on the cal- 
culated values of Ak s - The extinction map based on H-Ks colors 
is shown in the top left panel of Fig. [6] The top right panel of this 
figure shows the resolution of the extinction map, i.e. the radius 
at each position within which 20 stars were detected. The statis- 
tical uncertainty, RMS / V-V - 1 (with N = 20), of the extinction 
at each point in the map is shown in the bottom right panel of 
Fig. [6] It is < 0.1 mag for the entire map. As a cross-check, we 
also created an extinction map based on Ks-L' colors, in a com- 
pletely analogous way. It is presented in the bottom left panel of 
Fig. [6] and shows the same patterns as the extinction map based 
on H - Ks colors. Note that details will be more reliable in the 
map based on H - Ks because of the higher point source density 
measured in these bands. Based on the assumption that the ex- 
tinction law between Ks and L' is a power-law, we determined 
the corresponding power-law index, a^s-L', by requiring that the 
average Ak s must be the same, based on H - Ks or Ks - L' 
colors. A value of ans-n = 1-35 has been obtained, with a sta- 
tistical uncertainty of 0.12 and a systematic uncertainty - due to 
uncertainties in the effective wavelengths and overall photomet- 
ric calibration - of 0.29. This value of cuks-u agrees very well 
with the one derived from the location of the RC bump in the Ks 
and U luminosity functions. 

A histogram of the values of Ak s measured based on the 
H - Ks colors is shown in Fig. [7] The histogram is fit well by 
a Gaussian with a mean value of Ak s ,mean = 2.74 mag and a stan- 
dard deviation of 0.30 mag. 

The extinction map was used to produce a differential extinc- 
tion corrected color magnitude diagram of Ks vs. H - Ks. It is 
shown in Fig. [8] after exluding all stars with H-Ks < 1.8 as 
foreground and H-Ks > 2.8 as background stars. It can be seen 
how the extinction correction significantly reduces the scatter in 
the diagram (compare with the uncorrected diagram in the upper 
left panel of Fig.|4]i. The red clump is very well defined and the 
giant sequence can be perceived clearly to the right of the dashed 
line, which indicates the median H-Ks color. Stars on the left of 
the dashed line are expected to be either foreground or early-type 
stars. Known early-type stars (Buchho lz et al.|2009[ ) are marked 
with blue circles. They provide a sanity check on the reliability 
of identifying early-type stars via extinction corrected H - Ks 
colors. At Ks < 14 about 65% of the known early-type stars are 
identified correctly. How high is the probability to misidentify 
late-type stars as early-type with this method? Since the FOV of 
our work is larger than in comparable work, we limit the anal- 
ysis for investigating this question to the region within 12" of 
Sgr A* and to magnitudes brighter than Ks = 14. With these 
limits we can reasonably assume that the number of known and 
correctly classified stars (as early- or late-type) from literature is 
close to 100%. Within 12" of Sgr A* we find that about 30% of 
the stars brighter than Ks — 14 identified by the extinction cor- 
rected H-Ks colors may be erroneous identifications. Potential 
important sources of error are, for example, photometric errors 
and inaccuracies of the extinction map because of variability at 
scales below the resolution of our extinction map or because 
of locally different extinction related possibly with individual 
sources, e.g. due to circumstellar material. Nevertheless, consid- 



R. Schodel et al.: Peering through the veil: near-infrared photometry and extinction for the Galactic nuclear star cluster 



9 



ering the crudenes of the method presented here, a success rate 
of ~ 2/3 appears reasonable and suggests that our method may 
be useful for large-scale surveys in order to identify candidates 
for massive, young stars in highly extincted areas. At magnitudes 
Ks > 14 the extinction corrected H - Ks criterion is too crude 
and does not allow to separate early- and late-type stars in our 
sample. However, future work may provide more accurate colors 
and extinction values. 

We have marked all stars with Ks < 14 and a color bluer 
than the median H - Ks (after excluding foreground stars, see 
above) with crosses in the Ks image shown in Fi g,[Tj Many 
early- type stars reporte d by |Paumard et alT] ( |2006| l; |Buch holz 



et al. 



( 2009| l; Do et al.| ( |2009a) > and other authors are correctly 
identified. Following the error analysis in the preceding para- 
graph of the order 30% of the stars may be spurious identifica- 
tions and a similar number of early-type stars may be missed by 
this method. Nevertheless, this is the first time that early-type 
stars at the GC can be identified via broad-band imaging. Even 
taking into account the large estimated uncertainties, this rather 
crude method confirms the finding of Buchholz et al. ( 2009 ) that 
there exists a population of early-type stars beyond a projected 
radius of 0.5 pc from Sagittarius A*. 

Finally, we have used the extinction map to produce aKs lu- 
minosity function corrected for differential extinction (fore- and 
background stars excluded, mean extinction A Ks - 2.1 A, no in- 
completeness correction applied). It is shown in Fig. [9] By com- 
paring the power-law fit (dashed line) with the data, one can esti- 
mate that in the shallow observations used in this work (in order 
to avoid saturating the brightest stars) the completeness drops 
below 50% already at Ks « 17 (see also Table[2]l. 

6. Discussion 

For comparison of the derived photometry and extinction with 
other work, it is important to refer to an effective wavelength, 
i.e. to take into account such factors as the SED of the stars, ex- 
tinction, and the filter transmission curves. The measured value 
of the exponent of the power-law extinction-law is strongly de- 
pendent on the effective w avelength (see, e.g., |Nishiyama et aL 
|2006a|[Gosiing et al.|200"9"l|Stead & Ho are 2009). The transmis- 
sion curves of the NACO filters are available at the ESO website. 
In order to facilitate comparison of our work, we have calculated 
the effective wavelengths for the photometry presented here, us- 
ing equation (A3) of Tokuna ga & Vacca| ( |2005] ). 

We have found a significantly steeper extinction law at NIR 
wavelengths than in earlier work (e.g., Rieke & Lebofsky 1985 1 
Draine 1989; Mathis 1990, see the latter paper for a review ). 
There exists, however, very good agreement between the values 
for czh-ks = 2.21 + 0.24 (and an-Ks = 2.54 + 0.15, respectively) 
and a Ks _u = 1.34 + 0.29 (and a Ks _ L , = 1.35 + 0.31, respec- 
tively) found in this work and the respective values from other 
recent research. Gosling et al. (2009) report aj-Ks - 2.64 + 0.52 
based on NIR observations of the nuclear bulge. Note that they 
find variability of this power-law index along the line-of-sight, 
on scales as small as 5". We have not examined spatial variabil- 
ity of the power-law index in our work. Since we do not have 
/-band observations available and since the extinction law ap- 
pears to be clearly different between J and Ks and Ks and L', 
we cannot determine the extinction-law toward individual stars 



or small groups of stars like Gosling et al 
[eTaLl ( |2D59l i find a H - Ks = 2.0 and a Ks - 



(20091. Nishiyama 



3.6 



1.37 



based on 



NIR observations of the Galactic nuclear bulge. Stead & Hoare 
(2009) examined 8 regions from the UKIDDS Galactic Plane 
Survey and report q;nir = 2.14 ± 0.05. These values agree with 




Fig. 7. Histogram of the extinction, Ak s , based on H-Ks colors. 




Fig. 8. Differential extinction corrected color-magnitude dia- 
gram for the GC (assumed mean extinction is Ak s = 2.74). The 
dashed line indicates the mean H - Ks color. Known early-type 
stars are marked by blue circles. 



our results within the l<x - combined statistical and systematic - 
uncertainties. An exception is the value oh-ks - 2.54 + 0.15 de- 
termined from the color dependence of the RC in the H — Ks vs. 
Ks color-magnitude diagram. Nevertheless, this error lies within 
the 1 cr uncertainty of the value aH-Ks = 2.21+0.24 that we have 
derived in this work from the RC positions in the Ks and H lumi- 
nosity functions. We consider the latter value more trustworthy 
because it is more straightforward to determine (no cutting out 
of regions from the color-magnitude diagram), is based on large 
numbers of stars, and agrees very well with the recent works 
cited above. 

The extinction map based on H - Ks colors from this work 



resembles closely the one presented in Buchholz et al. (2009). 
The overall extinction we find in the present work is, how- 
ever, lower than the one reported by |Buchholz et al.| ( |2009| ) - 
A Ks = 2.1 A- ± 0.30 compared to A Ks = 3.14 + 0.43. The reason 
for this difference is the assumption of a flatter extinction law 
(based on |Draine|1989|) in|Buchholz et al.| ( |2009| ). The extinction 
map presented by Schod el et al.| ( |2007a[ ) is also similar to the one 
presented here, but shows a smaller FOV and is also based on a 
different extinction law. Agreement can also be seen between the 
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™<3 Ks 

Fig. 9. Ks luminosity function for the GC, fore- and back- 
ground stars excluded, and corrected for differential extinction. 
The mean extinction has been set to Ak s = 2.74. This luminosity 
function has not been corrected for completenes. The red dashed 
line is a power law fit to the stars 11 < Ks < 14. It has a power- 
law index of 0.27 + 0.02. 



extinction map produced in this work and the one presented by 
Scoville etdJpOOl) . The map in the latter paper encompasses 
a much larger region, but is less homogeneous than the one pro- 
duced here because it is derived from the line emission of ionized 
gas. Due to the strong radio emission from Sagittarius A* the 
latter authors could not reliably constrain the extinction in the 
immediate surroundings of Sagittarius A*. Scoville et al. ([2003 ) 
used an extinction-law of the form A A oc A~ L °, significantly dif- 
ferent from what is found here and in other recent work on this 
subject. 

A caveat may apply to the regions with the strongest extinc- 
tion. As can be seen in Fig.[T[ there are some regions (e.g. at an 
offset of ~ +12" E and ~ -13" S of Sagittarius A*) with very 
high extinction and consequently few stars. Extinction may be 
under-estimated by our method in these patches of extremely 
high extinction because the number of stars that can be detected 
behind a strong extinction screen may be very small. 

Finally, useful for analyses of the so-called "S-cluster" of 
stars that are tightly bound to Sagittarius A* and of the NIR 
emission of SgrA* itself, we can use the extinction map pre- 
sented in Fig. [6] to calculate the mean extinction and its stan- 
dard deviation toward Sagittarius A*. For a circular region with 
a radius of 0.5" centered on Sagittarius A* we obtain A^sgrA* = 
4.35 ± 0.12, A Ks , SgrA , =2.46 + 0.03, andA^, Sg , AH , = 1.23 + 0.08. 
The uncertainty of A Ks y s gr A* corresponds to the standard devia- 
tion of extinction within the circular aperture. The uncertainties 
for A[],sgrA* and Au t s gr A* include additionally the uncertainties 
of (Xh-ks an d &Ks-L'- These extinction values are different from 
the ones that have typically been used in the analysis of NIR 
lightcurves from Sgr A*, e.g., A H = 4.3 , A K = 2.8, A v = 1.8 
(|Genzel et al.| [2003 ; Ec kart et al. 
2009b|), or A L ' = 1.8 ( |Ghez et aL 



2006), A K = 3.3 flDo et al.| 



2004| >. While the impact on 



emission/accretion models of SgrA* will probably be minor, we 
point out that the new extinction values at Ks and L' imply that 
the de-reddened emission from Sgr A* may be weaker than as- 
sumed by a factor 1.4 - 2.2 at Ks and a factor of ~ 1.7 at L' . 



A correct estimate of the extinction is crucial to address one 
of the long standing issues in the central parsec, namely the ap- 
parent incompatibility between the ionizing flux arising from the 
massive stars and the ionisation of the gas (Najarro et al. 1997; 



Lutz 1999[ ). Further, these authors found that the population 



of Hel stars, displaying fairly large luminosities was unusually 
high. Such incompatibility was recently claimed to be solved by 



Martins et al. ([2007 ) with the inclusion of line-blanketing in the 



models, leading to lower bolometric luminosities and lower UV 
fluxes for the Hel stars. However, we attribute most of the dis- 
crepancy between the results of Najarro et al. ( 1997 ) and M artins 
et al.| ( |2007l l for the Hel stars to the estimate of Ak between both 
investigations. Indeed, the final Ak values computed by |Najarro 
1997 i were 50% larger than the ones adopted by 



etal. 



et al. 



Martins 



(2007) from the extinction map of Schodel et al. ( 2007a I. 



Such difference was translated into roughly one magnitude in lu- 
minosity and was responsible for the much lower luminosity of 



the Hel stars revised by Martins et al. (2007 1. Similarly, the na- 
ture of the "S" stars is also affected. While Marti ns et"aT1 ([2008) 
adopted Ak = 2.25 for S2 it turns out that the extinction map 
they used was calibrated on S2 with an Ak — 2.8 Schodel et al. 
( 2007a| l. That would have changed S2 from being a dwarf to be- 
ing a giant. It turns out that the actual value of Ak s = 2.46 toward 
S2 is not too different from the value assumed (erroneously) by 
Mar tins et al.| (2008 ). As we show in this work, and as other sim- 
ilar recent investigation has shown, the solution and key to these 
problems is that the power-law slope of the NIR extinction law 
is significantly steeper than what was assumed in earlier work. 
This leads to lower extinction in the /f-band, but significantly 
stronger extinction at shorter wavelengths. 




Fig. 10. Ilustration of the relation between the stars detected at 
Ks (blue), the Ak s extinction (green), and the diffuse emission 
in the mid-infrared (PAH I, red). 



In Fig.[l0]we show a color composite image of the stars, de- 
tected at Ks, the A Ks extinction, and the MIR (PAH I) emission 
(data from |Schodel et al.|2007b) . It can be seen how the extinc- 
tion is low in a region running roughly from SW to NE. This ere- 
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ates an apparent high stellar density and therefore elongation of 
the NSC along the galactic plane (as was also shown by Schodel 
et al. 2007a|l. Similarly, the lower extinction to the NW of Sgr A* 



is responsible for the higher stellar density observed in that re- 
gion. The strong influence that extinction has on the appearance 
of the nuclear star cluster can also be nicely seen in Fig. 1 of 
|Scoville et al.| ( (2003] >. There appears to be no obvious relation 
between the warm dust emitting in the mid-infrared and the ex- 
tinction. This is no surprise because the surface density of HI and 
H2 and the warm dust associated with it in the mini-spiral is just 
of the order a few xlO 4 cm" 3 (see, e.g., discussion of this topic in 
Moul taka et al.||2009] ). The extinction toward the GC is largely 
produced by clouds in the foreground, with the possible excep- 



tion of some dust-enshrouded individual objects (Moultaka et al. 
[2009] l . 

For our analysis we have used a weighted average of recent 
measurements of the distance to the GC, assuming that all those 
measurements are independent. We did not investigate this as- 
sumption in depth. However, we took into account the system- 
atic uncertainties of the individual measurements when comput- 
ing the weighted average. As a further simple test we can use 
sub-samples of the measurements and check the resulting values 



of Rq for consistency. Using only geometrical estimates ( Ghez 



et al.|2008HTrippe et al.|2008"l|Gillessen et al.|2009[|Reid et al 



2009a b) we obtain Ro, ge om,i - 8.20 + 0.22 kpc (emitting the esti- 
mate from the S0-2/S2 orbit by |Ghez et al.|2 008 ) and Ro, ge0 m.2 = 
8.10 + 0.26 kpc (emitting the estimate from the S0-2/S2 orbit by 
Gille ssen et al.|2009] l. Using only measurements based on vari- 
able and RC stars, we obtain Ro, s i„ rs = 7.85 ± 0.24 kpc. As an 
additional test we computed the unweighted average of all mea- 
surements: Ro,unweighted - 8.05 + 0.28 kpc. It agrees within the 
uncertainties with the weighted average. If we remove the mea- 
surement with the strongest apparent deviation (7.52 ± 0.36 kpc, 
|Nishiyama et al.|200 6b ), we obtain Ro tm weighted - 8. 12+0.20 kpc. 
The unweighted averages of the sub-sets of distance measure- 
ments considered above agree also well with the weighted re- 
sults. We therefore believe that the value for/?o = 8.03+0.15 kpc 
derived in this work and its related uncertainty can be considered 
reliable. 

We believe that future work should focus on deriving the 
power-law index for the NIR extinction-law toward the GC with 
higher precision and examine its variability over the central par- 
sees. This could be achieved, for example, by deep AO imag- 
ing using a large number of narrow band filters across the near- 
infrared. The main limitation for this kind of work may come 
from the still limited performance of AO systems at short NIR 
wavelengths (/-band), combined with the long exposure times 
required by the high extinction. An additional great difficulty, in 
fact the main source of uncertainty in this work, is the difficulty 
to obtain precise and accurately calibrated photometry from AO 
data (crowding, large seeing foot of PSF, spatial variability of 
PSF). We have, however, no doubt that there will be significant 
technical progress in AO techniques and related data processing 
software in the next few years. 



7. Summary 

We provide photometry for ~ 7700 point sources extracted from 
adaptive optics observations of the central parsec of the Galactic 
center in the Ks, and for many of them also in the H and L' 
bands. The data are made publicly available as electronic on- 
line material. Adaptive optics photometry has been published 
previously in some papers. However, we believe that the data 
presented in this work supersede those earlier results or are 



complementary to them because of reasons such as: (a) limited 
field-of-view of previous observations or focus on the central 
few 0.1 pc; (b) isoplanatic effects were not taken properly into 
account; (c) bias by the application of methods such as Lucy- 
Richardson deconvolution (the companion paper to this work, 
by Schoedel, A&A, 2009, accepted for publication, explains our 
photometric method and points out the bias introduced by the 
Lucy-Richardson method); (d) no transparent explanation of the 
photometric calibration. This data set is the most comprehensive 
one published so far. 

Based on the photometric results, we analyze extinction to- 
ward the central parsec of the GC. Absolute values of extinction 
in the three examined bands are derived based on the well-known 
properties of red clump stars, without the need to measure stellar 
colors and without having to assume an extinction law. The val- 
ues obtained are A H = 4.48 + 0.13 mag, A Ks = 2.54 ± 0.12 mag, 
and Au = 1.27 ± 0.18mag. This results in Ah : Ak s ■ An = 
1.76 : 1 : 0.50. Assuming the validity of a power-law ex- 
tinction law between H and Ks and between Ks and L', this 
is equivalent to power-law indices of an-Ks - 2.21 + 0.24 
and aKs-L' = 1.34 ± 0.29. This is significantly steeper than in 
earlier work (e.g. Rieke & Lebofsky 1985) |Draine| fl989), but 
agrees with recent studies of extinction toward the nuclear bulge 



( Nishiyama et al.| 2009 Gosling et al. 2009 1 and other targets 
in the Galactic plane ( [Stead & Hoare 2009). We also confirm 
the flattening of the extinction law beyond ~ 3fim (see also 
I Lutz et al.|T9"96l|Indebetouw et al.|2005i|Viehmann et al.|2003j 
Flaherty et al.|2007||Nishiyama et al.|2009| >. 

Based on H—Ks colors of the stars, we provide a detailed ex- 
tinction map for a FOV of 40" x 40" (~ 1 .5 pc x 1 .5 pc) centered 
roughly on Sagittarius A*. The spatial resolution of the map is ~ 
1". Both statistical and systematic uncertainties of the extinction 
map are smaller than 10%. The extinction map is made publi- 
caly available as electronic on-line material. We find that extinc- 
tion varies significantly on arcsecond scales. The mean extinc- 
tion based on stellar colors is found to be Ak s = 2.74 + 0.30 mag, 
in good agreement with the value of mean extinction based on 
the RC method. An extinction map based on Ks - L' colors, al- 
though less accurate because of lower point source surface den- 
sity and possible contamination by diffuse emission from dust in 
the mini-spiral, confirms the patterns seen in the extinction map 
based on H — Ks colors. Requiring the mean extinction between 
the two maps to be the same, we derive (Xks-l 1 = 1.35 + 0.31, 
in excellent agreement with the value of this parameter derived 
from the mean magnitudes of the RC bump (see above). 

Applying the derived extinction map to the observed H — Ks 
colors of the stars, we are able to correct the corresponding 
color-magnitude diagram. This reduces the scatter in the dia- 
gram considerably and makes it possible to distinguish between 
stars in the giant branch and young, massive stars. We confirm 
recent findings of the presence of of young, massive stars at pro- 
jected distances > 0.5 pc from Sagittarius A*. 

A minor, but very interesting by-product of this work is a 
weighted average of recent (years 2006-2009) estimates of Rq. 
We find that the distance of the GC is constrained with an uncer- 
tainty of just 2%: R = 8.03 ± 0.15 kpc. 

For a circular region with a radius of 0.5" centered on 
Sagittarius A* we obtain A HySg ,A* - 4.35 ± 0.12, A K ^ SgrAt - 
2.46 + 0.03, and A L , >Sgr A* = 1.23 + 0.08. These values imply that 
the NIR emission from Sgr A* in the Ks and L' -bands may be 
weaker than previously assumed, by factors 1 .4 - 2.2. 

We have used the following effective wavelengths in this 
work: A H = 1.677 ± 0.018/mi, A Ks = 2.168 + 0.012/mi, and 
^e«,u - 3.636 + 0.012yum. Note that the effective wavelength 
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is a complex function of the transmission functions of the at 
mosphere, the filter, the optical system, extinction, and stellar 
type. For discussions of this problem, see, for example, Espinoza 
|et al.| ( |2009) and |Stead & Hoare| ( |2009] l. The uncertainties cited 
for the effective wavelengths are based on best-estimates of the 
variablity of these values for the case of Galactic center observa- 
tions and the observing conditions and instruments used in this 
work. 
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Appendix A: Effective wavelength 

We calculated the effective wavelengths of our observations fol- 
lowing equation (A3) of Tokuna ga & Vacca] ( |2005| l. The trans- 
mission curves for the NACO H and Ks filters were downloaded 
from the instrument web sit^] Spectra of the atmospheric trans- 
mission were downloaded from the Gemini telescope web sitq^ 



They are based on models with the ATRAN software (Lore 
1992). The transmission curves for 3 mm precipi table water va- 
por were used because humidity was high during the H and Ks 
observations. 

We did not use any particular stellar atmosphere model for 
the calculations. Instead, simple extincted blackbody models 
were used. Table|A| lists the values of the effective wavelength 
for various values of Ak and effecti ve temperatures of the b lack- 
bodies. A power-law Ak °c At" (Nishiyama et al. 20Q9| > was 
used to calculate the extinction at different wavelengths. We as- 
sumed ctH-Ks - 2.22 for the extinction between H and Ks , and 
®Ks-L! - 1.33 for the extinction between Ks and L' ^\ The ma- 
jority (> 90%) of stars in the images analyzed in this work are of 
late-type, with a large percentage of RC stars (see, e.g., Buchho 
et al.|2009] ). We therefore adopt as effective wavelengths /l e ff,H = 
1 .67/ /mi, /tefl^Ks = 2.168 /mi, and /i e ff,L' = 3.636/mi for an aver- 
age extinction of A ks ~ 2.5. 

The uncertainty introduced by the choice for the amount 
of precipitable water vapor is low. Choosing 1 .6 mm instead of 
3.0mm will change the effective wavelengths by < 0.003 //m 
for H, < 0.001 /mi for Ks, and < 0.01 /mi for V . Changing 
the power-law index for the extinction curve between H and Ks 
to an-Ks — 2.0 (Nishiyama et al. 2009) results in a change of 

< 0.003 //m for /l e ff,Ks, and < 0.009 /mi for A e g tli . Changing the 
power-law index for the extinction curve between Ks and L' to 
&Ks-L' = 1-37 (Nishiyama et al. 2009| l results in a change of 

< 0.001 /mi for /t e ff,L'- The uncertainty in extinction introduces 
an uncertainty of < 0.015 /mi for H, < 0.01 /mi for Ks, and 

< 0.002 /mi for V . The unknown effective temperature of the 
stars introduces uncertainties of < 0.005 //m for H, < 0.006 //m 
for Ks, and < 0.001 /mi for U. 

We estimate the combined 1 cr uncertainty of the effective 
wavelengths to 0.018 //m for H, 0.012 /zm for Ks, and 0.012 /mi 
fori'. 



3 http://www.eso.org/sci/facilities/paranal/instruments/naco/index.html 

4 http ://w w w. gemini . edu/node/ 1 07 8 1 ?q=node/ 10789 

5 This value of ctKs-L' is from a preliminary estimate at the beginning 
of this work and is smaller than the value a Ks -u = 1.38 that we derived 
in the end. However, the associated change of /t e tf,L' is minimal, see 
discussion of the uncertainties in the next paragraph. 
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Table A.l. Effective wavelengths in units of //m for the NACO 
H, Ks, and L' filters for blackbodies with different effective tem- 
peratures and extinction, A Ks . 



H 






2.0 a 2.5 3.0 3.5 


3000" 
4700 
30000 


1.674 1.681 1.688 1.695 
1.669 1.677 1.684 1.691 
1.664 1.672 1.680 1.687 


Ks 




3000 
4700 
30000 


2.167 2.171 2.176 2.180 
2.163 2.168 2.173 2.177 
2.159 2.164 2.169 2.174 


L' 




3000 
4700 
30000 


3.636 3.637 3.638 3.639 
3.635 3.636 3.637 3.638 
3.634 3.635 3.636 3.637 



a Extinction in the Ks band, A Ks . 
b Effective temperature. 
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Table A.2. List of detected point sources in the NACO H, Ks, and L'-band observations. We only include the first ten lines of the table in the 
printed edition of this work. 



ID 


R-SgrAi.'' 
["] 


R.A. b 
["] 


Dec. c 
["] 


t rd 

H a 


AH e 


Ks 1 


aa:.v e 


L " 


AL" 


1 


0.028 


0.009 


-0.027 


17.50 


0.05 


15.7 


0.04 


0.0" 


0.0 


2 


0.107 


0.031 


-0.102 


18.19 


0.06 


16.39 


0.04 








3 


0.163 


-0.006 


0.163 


16.00 


0.02 


14.13 


0.01 


12.84 


0.03 


4 


0.231 


0.189 


0.132 


16.93 


0.03 


15 


0.01 








5 


0.245 


-0.231 


-0.081 


18.69 


0.07 


16.75 


0.05 








6 


0.265 


0.022 


-0.264 


16.53 


0.03 


14.7 


0.02 








7 


0.321 


0.296 


0.124 


16.47 


0.02 


14.61 


0.01 








8 


0.347 


-0.274 


0.213 


17.45 


0.03 


15.66 


0.09 


14.35 


0.06 


9 


0.354 


-0.053 


0.35 


17.22 


0.03 


15.39 


0.02 








10 


0.363 


0.02 


-0.363 


16.23 


0.03 


14.12 


0.01 


12.6 


0.03 



Distance from Sgr A* projected on the sky. 

Offset in right ascension from Sgr A*. Uncertainties in the astrometric positions can be as large as ~ 0.1" because high precision astrometry 
was not at the focus of this work. The conversion to offsets from Sgr A* was done by assuming a pixel scale of 0.027" per pixel (ESO NACO 
manual). The camera rotation angle was assumed Odeg east of north and no camera distortion solution was applied. 
Offset in declination from Sgr A*. For more information see also footnote b . 
Magnitude in the NACO H band filter. 

Uncertainty in H, obtained by quadratically combining the formal uncertainty from PSF fitting, based on Gaussian read-out noise and Poisson 
photon noise, and the uncertainty due to the limited knowledge of the PSF. The 1 <x uncertainty of the zero point, AZP H = 0.06, is not included 
in this column and has to be taken into account when considering absolute photometry. 
Magnitude in the NACO Ks band filter. 

Uncertainty in Ks. For more information see also footnote c . The 1 tr uncertainty of the zero point is AZP Ks = 0.06. 
Magnitude in the NACO L' band filter. 

Uncertainty in L' . For more information see also footnote e . The 1 <x uncertainty of the zero point is AZP L > =0.15. 
A value of 0.0 in the table indicates that the corresponding measurement is not available. 



